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ABSTRACT 


This thesis presents two methods for determining the bearing of a source 
generating a Very High Frequency (VHF) bridge-to-bridge radio transmission. Using 
principles of interferometry, one can utilize the calculated phase difference between the 
received signal at multiple antennas to determine the Angle of Arrival of the detected 
transmission. This translation of phase difference to Angle of Arrival is accomplished 
using equations based on signal properties and geometry. The theoretical method is 
shown proving the relationship between Angle of Arrival and phase difference, as well as 
how a single platform could accomplish this detection and calculation. Theoretical 
simulation was accomplished using various simulation tools including Mathwork’s 
Simulink and Tonne Software’s Elsie. Methods are then provided to detect the phase 
difference using both a series of analog mixers and filters as well as digitally, using 
software radio. Analog filters were built and tested to determine the relationship between 
phase difference and voltage output. Software programs were written for a Software- 
Defined Radio implementing digital filtering to verify the analog performance. Results 


and accuracy are shown based on initial testing. 
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EXECUTIVE SUMMARY 


Situational awareness at sea is essential for the navigational safety and security of 
the United States’ military vessels. The Navy routinely operates in high-traffic areas with 
large volumes of shipping, including straits, canals, and other waters close to shore. As 
the volume and shipping density increases, increasing coordination is required between 
vessels, both military and civilian, to ensure the safety of all. Much of that coordination 
is conducted over Very High Frequency (VHF) bridge-to-bridge Radio. This maritime 
system operates between 156 and 162 MHz with 50 kHz of bandwidth assigned to each 
channel. It is by using this system that ships hail one another and conduct conversations 
regarding maneuvering intentions and other safety issues. In areas of high traffic, there 
can be nearly continuous chatter on this system. Unless one is participating in the 
conversation, it can be very difficult for many ships to determine who is talking at the 
time and from where the conversation is originating. Even if not involved in the 
conversation, it could be vital that a United States naval vessel be aware of the 
maneuvering intentions of the ships around it, for its own safety and security. Most 
conversations on bridge-to-bridge radio either originate or are conducted on channel 16, 
which is located at 156.8 MHz. This channel is reserved for international distress, safety, 
and calling; most ships and coastal stations maintain a listening watch on this channel. 
This thesis focuses on this frequency, but the methods conducted could be used to 


achieve similar results on any channel in the bridge-to-bridge VHF frequency range. 


Previous research into the area of direction finding has focused on using 
techniques such as Time Difference of Arrival (TDOA) and Frequency Difference of 
Arrival (FDOA). The drawback to these systems is that they either require multiple 
transmitting platforms or multiple receiving platforms. While United States naval vessels 
traditionally operate in groups, there are situations where they are transiting alone. In 
these situations, many ships would benefit from a direction-finding system that can 
obtain a solution with only one transmitting platform and one receiving platform. 
Furthermore, the research focuses on simpler, low cost solutions to this problem. This 


can be accomplished using multiple antennas on a single platform and measuring the 


Xlil 


phase differences seen at each antenna. This technique will enable a single ship to 


determine the bearing of vessels transmitting on their bridge-to-bridge radios. 


This thesis examined two possible solutions to this problem, an analog solution 
and a digital solution. The analog solution was explored using traditional Radio 
Frequency (RF) components: resistors, capacitors, and inductors. These were utilized 
along with traditional RF mixers to achieve a solution. Rapid improvements in 
technology, however, have led to many functions traditionally found in hardware to now 
be implemented in software. In many cases, software-based solutions can be preferable 
as they can be cheaper, more flexible, and achieve higher performance. It is these 
advantages that led to the development of the software radio. There are exceptions, of 
course; in this research, the hardware-based solution proved cheaper than the software 


solution. 


Ideally, the software-defined radio receiver would be nothing more than an 
analog-to-digital converter (ADC) attached to an antenna. The ADC would digitize the 
analog signal, and then computer software would process all the results and perform any 
needed applications or calculations. However, due to the desired frequencies of the 
signals of interest, oftentimes a downconverter is required to bring the signal to a lower 
intermediate frequency or baseband before digital sampling can occur. The device 
utilized in this thesis to do both downconversion and analog to digital conversion is the 


Universal Software Radio Peripheral (USRP). 
The objectives of this research were to: 


e Review existing methods of using phase difference of a received signal as 


an inexpensive means of determining a line of bearing to the transmitter. 
e Model the overall system using computer simulations. 
e Construct circuit boards based on design results 
e Test the analog system to verify performance 


e Use software to achieve a digital solution 
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e Test the digital system to verify performance 


e Compare analog and digital results to determine the preferred solution to 


the research problem 


These objectives of this research were met with different levels of success for the 


analog and digital systems. 


The feasibility of using phase difference as a means of obtaining the Angle of 
Arrival was verified by mathematical calculations, which showed a direct link between 
the phase difference of one signal received at two antennas and the Angle of Arrival from 
which the signal originated. The system was then implemented in Simulink® to 
determine the feasibility of implementing the mathematical solution. While the Simulink 
simulation proved highly successful, its direct implementation in hardware proved 
infeasible due to the high frequency signal (156.8 MHz) and the relatively narrow 
bandwidth (50 kHz). The design of the system was then altered to take into account the 
hardware limitations, and the system was re-simulated. The system proved successful, 


which led to the design of the hardware. 


The design required a series of analog filters which were then modeled in Elsie to 
ensure feasibility with industry’s currently manufactured RF component ratings. Once 
the design was complete and verified, circuit boards were constructed using commercial 
off-the-shelf RF components and surface-mounted soldering to circuit boards. These 
boards were tested successfully and demonstrated a clear relationship between the 


voltage output of the system and the phase difference. 


Concurrently, a digital system was designed using the GNU radio and Octave 
software projects. In this system, all filtering was performed digitally and the output was 
sent to a personal computer where the results were displayed. As in the analog case, 


there was a clear relationship seen between the output and the phase difference. 


To further make this system compatible for shipboard use, the following 


additional research is required: 


XV 


e Test using wireless signals to determine system effectiveness. 


e Implement Low Noise Amplifiers and Automatic Gain Control to 


account for real-world environmental conditions. 
e Test in a shipboard environment to determine maritime suitability. 


With free and open sea lines of communication remaining vital to the United 
States economic and vital military interests, further research is recommended, for this 
research shows excellent potential for an inexpensive solution to the bridge-to-bridge 
radio direction-finding problem, thereby increasing situational awareness for United 


States warships. 
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I. INTRODUCTION 


A. BACKGROUND 


Situational awareness at sea is essential for the navigational safety and security of 
the United States’ military vessels. The Navy routinely operates in high-traffic areas with 
large volumes of shipping, including straits, canals, and other waters close to shore. As 
the volume and shipping density increases, more coordination is required between 
vessels, both military and civilian, to ensure the safety of all. Much of that coordination 
is conducted over Very High Frequency (VHF) bridge-to-bridge Radio. This maritime 
system operates between 156 and 162 MHz with 50 kHz of bandwidth assigned to each 
channel [1]. It is by using this system that ships hail one another and conduct 
conversations regarding maneuvering intentions and other safety issues. In areas of high 
traffic, there can be nearly continuous chatter on this system. Unless one is participating 
in the conversation, it can be very difficult for many ships to determine who is talking at 
the time and from where the conversation is originating. Even if not involved in the 
conversation, it could be vital that a United States naval vessel be aware of the 
maneuvering intentions of the ships around it, for its own safety and security. Most 
conversations on bridge-to-bridge radio either originate or are conducted on channel 16, 
which is located at 156.8 MHz. This channel is reserved for international distress, safety, 
and calling; most ships and coastal stations maintain a listening watch on this channel. 
This thesis focused on this frequency, but the methods conducted could be used to 


achieve similar results and other channel in the bridge-to-bridge WHF frequency range. 


Navigational safety was not the only relevant concern prior to commencing work 
on this thesis. A proposed solution for using bridge-to-bridge radio as a direction-finding 
tool could help provide added safety and security to many ships and their crews. A major 
motivator to this research was an incident that occurred in the Strait of Hormuz on 
January 6, 2008. On that day, three United States Navy ships were completing a routine 
transit through the strait when they were aggressively approached by five Iranian small 
boats. It was over bridge-to-bridge radio that the small boats communicated their threats 


against the U.S. ships, telling them that they would soon explode [2]. In this scenario, the 
1 


warships did have visual contact with the aggressive small boats, but a situation could be 
envisioned where threats are received from ships outside visual range or are received 
from an area with high traffic density. In littoral environments, the signals could also be 
originating from shore-based installations. The only way to tell who is transmitting is to 
trust vessels to identify themselves truthfully or to have a system to map signals to their 
sources. In these situations, a system utilizing bridge-to-bridge radio as a direction- 
finding tool would be of great benefit to many warships. This system would give a line 
of bearing in the direction from which the threats are being transmitted. In a situation 
where no other ships are in sight, it would give the warship extra time to get ready for the 
threat and a direction as to where it should focus its efforts. In a heavy shipping area, it 
would help the warship distinguish potential threats from the mass of surrounding traffic. 
In both cases, the safety and security of the warship would be enhanced, making it better 


prepared for any potential threat. 


This thesis examined two possible solutions to this problem, an analog solution 
and a digital solution. The analog solution was explored using traditional Radio 
Frequency (RF) components: resistors, capacitors, and inductors. These components 
were utilized to create filters that, combined with traditional RF mixers, achieved a 
solution. Rapid improvements in technology, however, have led to many functions 
traditionally found in hardware to now be implemented in software. In many cases, 
software-based solutions are preferable as they can be cheaper, more flexible, and 
achieve higher performance. It is these advantages that led to development of the 
software radio [3]. The primary disadvantage to a software solution is the relatively 
recent development of this field. It may take time to develop a software solution to a 
problem that has traditionally been solved with hardware. The final solution with 
software may be preferable when compared with hardware, but the time to develop that 


solution may prove problematic if a strict timeline exists. 


Ideally, the software-defined radio receiver would be nothing more than an 
analog-to-digital converter (ADC) attached to an antenna. The ADC would digitize the 
analog signal, and then computer software would process all the results and perform any 
needed applications or calculations. However, due to the desired frequencies of the 
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signals of interest, oftentimes a downconverter is required to bring the signal to a lower 
intermediate frequency or baseband before digital sampling can occur. The device 
utilized in this thesis to do both the downconversion and analog to digital conversion is 


the Universal Software Radio Peripheral (USRP) [3]. 


This research explored whether a single platform would be able to perform 
direction finding via an inexpensive and simple solution and whether the phase difference 
that occurs between one signal arriving at two antennas could be utilized for that purpose. 
This thesis sought to answer both those questions through both theory and 
experimentation. 


B. OBJECTIVES AND APPROACH 


The objectives of this research were to: 


e Review existing methods of using phase difference of a received signal as 


a means of determining a line of bearing to the transmitter. 
e Model the overall system using Simulink 


e Model the needed analog filters using Elsie, Tonne Software’s electrical 


filter design and analysis program 
e Construct circuit boards based on design results 
e Test the analog system to verify performance 
e Use GNU Radio and Octave programs to achieve digital solution 
e Test the digital system to verify performance 


e Compare analog and digital results to determine the preferred solution to 


the research problem 


The feasibility of using phase difference as a means of obtaining the Angle of 
Arrival was verified by mathematical calculations, which showed a direct link between 
the phase difference of one signal received at two points and the Angle of Arrival from 


which the signal originated. The system was then implemented in Simulink to determine 


the feasibility of implementing the mathematical solution. While the Simulink simulation 
proved highly successful, its implementation in hardware proved infeasible due to the 
high frequency signal (156.8 MHz) and the relatively narrow bandwidth (50 kHz). The 
design of the system was then altered to take into account the hardware limitations, and 
the system was re-simulated. The system proved successful, which led to the design of 


the hardware. 


The design required a series of analog filters which were then modeled in Elsie to 
ensure feasibility with industry’s currently manufactured RF component ratings. Once 
the design was complete and verified, circuit boards were constructed using commercial 
off-the-shelf RF components and surface-mounted soldering to circuit boards. These 
boards were tested successfully and demonstrated a clear relationship between the 


voltage output of the system and the phase difference. 


Concurrently, a digital system was designed using GNU radio and Octave 
programming [3, 4. In this system, all filtering was performed digitally and the output 
was sent to a personal computer where the results were displayed. As in the analog case, 


there was a clear relationship seen between the output and the phase difference. 


C. RELATED WORK 


The issue of direction finding as a research field is not new. Much research has 
been previously conducted exploring various techniques and methods to find a working 
solution. In 2006, Derek Elsanesser published a paper detailing his Discrete Probability 
Density Method for geolocation [5]. His solution was well suited for geolocation 
utilizing large numbers of measurements from different locations, for it effectively fused 
the data from multiple sensors. In this research, however, the ‘sensor’ utilized is the 
already-existing shipboard bridge-to-bridge radio receiver. Additionally, the 
measurements are not received from multiple platforms, but rather from two antennas on 


one platform. 


Other research into direction finding focused on using information other than 
phase to achieve the geolocation solution. Three popular techniques are: Angle of 
Arrival (AOA), Frequency Difference of Arrival (FDOA), and Time Difference of 
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Arrival (TDOA). Geolocation using Angle of Arrival requires lines of bearing taken at 
multiple locations at a single instance of time. Triangulation is then used to generate a 
solution from the individual lines of bearing. Time Difference of Arrival utilizes cross- 
correlation or some other delay-estimation technique from two spatially separated, 
moving receivers to generate a geolocation solution. Frequency Difference of Arrival is 
similar to TDOA but utilizes the difference in Doppler shifts between two spatially 
separated receivers to achieve the geolocation solution. There has been significant 
research into different ways of fusing these methods together to achieve a more precise 


geolocation solution [6, 7, 8]. 


The use of phase difference as a source of direction-finding is known as 
interferometry. It has been used as a common technique in electronic warfare direction- 
finding theory [9]. This thesis was concerned with the scenario of a single vessel 
steaming alone and needing to determine the location of a transmitter over bridge-to- 
bridge radio. This solution would also prove to be very cost-effective as the new 
equipment required would be minimal for it uses equipment already installed on the vast 


majority of marine vessels. 


Recently, many students at the Naval Postgraduate School have written theses on 
software-defined radio (SDR), but only a few have explored SDR as a method for 
direction finding. One thesis that proved particularly useful was Ian Larsen’s 2007 thesis 
which used SDR for direction finding [10]. However, his research focused on correlation- 
based Time Difference of Arrival to achieve a solution. As stated previously, this thesis 
focused on both hardware and software solutions with the signal’s phase difference as the 


basis of the solution. 


D. THESIS ORGANIZATION 


The body of this thesis is divided into three chapters to indicate the level of 
technical detail. Chapter II details the analysis and derivation of the phase difference- 
based solution to the direction-finding problem. Chapter III focuses on the hardware 
solution to this problem, including the design, experimentation, and results of the analog 


hardware solution. Chapter IV details the software solution to this problem, including 
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general information on Software-Defined Radio, the Universal Software Radio Peripheral 
(USRP), and the results obtained. Chapter V summarizes results, makes conclusions, and 


recommends future work. 


I. ANALYSIS 


This chapter describes the theory of interferometry and its mathematical basis. 
Mathematical relationships are then developed which show the basis for calculating the 


phase difference, which will be used to calculate the Angle of Arrival. 


A. INTERFEROMETRY 


As stated earlier, the basis for this research’s direction-finding solution is the 
phase difference that exists due to the arrival of a signal at two different points. This 
method of using the difference in phase to determine the Angle of Arrival of an incoming 
signal is known as interferometry. Figure 1 shows a graphical representation of the 


interferometry solution. 





Figure 1. Graphical Representation of Interferometry. (From [9]). 


The direction to the transmitter is labeled as AOA. The distance between the two 
receiving antennas (baseline) is defined as B. d is the distance between an antenna and 


the incoming wavefront. The mathematical relationships between these variables is [9] 


AOA = arcsin(4/,). (2.1) 


To calculate d, one uses the formula [9] 


__ ¢ 
er (2.2) 


8 
where c is defined as the speed of light (= 3.00 x10 n), f is the frequency of the 


transmitted signal, and @ is the phase difference (in degrees) of the received signal 


arriving at the two antennas. 


Using these formulas, one can see that the more phase change required to change 
the Angle of Arrival, the more accurate the direction-finding system will be. In other 
words, greatest accuracy will be achieved for angles near perpendicular to baseline and 
worst at angles near ends of the baseline. Additionally, accuracy will improve as the 
length of the baseline, relative to the transmitted signal’s wavelength, is increased. This 


is true up to a baseline equal to one-half of the wavelength. After that point, the phase 


difference changes more than 360° as the AOA moves from +90° to —90° [9]. 
Rearranging Equation (2.1) yields 
d = Bsin(AOQA). (2.3) 


Combining (2.2) and (2.3) yields 


Bsin(AQA) -—— 


AOA = arcsin Lie . 
360-B: f. 


(2.4) 


The above equation proves that the Angle of Arrival of the transmitted bridge-to-bridge 
signal can be calculated if one knows the phase difference encountered by the two 
separate antennas. The value of c is a known constant, and the variables, B and f, are 
known, for the user will know the distance between the receiving antennas, as well as the 
frequency of the signal one is capturing. Once the phase difference is known, one must 
use that information to calculate the Angle of Arrival. Angle ambiguities arise in that 


calculation, which will be detailed further in Chapter III.B.1.c. 
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B. CALCULATING PHASE DIFFERENCE 


As stated earlier and seen in (2.4), the key to determining the Angle of Arrival is 
knowing the phase difference of the incoming signal. This can be accomplished by 
multiplying the incoming signal with its phase-shifted version and then performing low- 
pass filtering. With the use of Automatic Gain Control and Low Noise Amplifiers, one is 
able to set the amplitude of the signals seen by the system. A cosine wave represents the 
carrier wave of the bridge-to-bridge signal. Assuming unit amplitude, the incoming 


signal is given by 
cos(2z f.t) (2.5) 
and the phase shifted signal is given by 
cos(2af.t+¢). (2.6) 
Multiplying the two signals together, with output a(t), yields 


a(t) =| cos(27 f+) || cos(27 fr) | 


(2.7) 
a(t) = 5008 pt 5008 (4c f.t+¢) 


Applying a Low Pass Filter to remove the high-frequency components yields the 


following relationship: 
1 
LPF {a(t)} = 3008 p (2.8) 


This research focused on incoming signals from a bridge-to-bridge radio system. 
Because a ship is a relatively slow-moving vessel, the phase difference will not be 
changing rapidly. Due to this fact, the output value of (2.8) is essentially a DC value. 
Thus, that value can be applied to (2.4) to calculate the Angle of Arrival. As stated 


earlier, angle ambiguities arise when @ is used to calculate Angle of Arrival, which will 


be detailed in Chapter II.B.1.c. 


C. CALCULATING PHASE DIFFERENCE USING LOCAL OSCILLATOR 


In the hardware solution, which will be detailed further in Chapter III, it proved 
necessary to introduce a local oscillator to calculate the phase difference at the two 


receiving antennas. The system with the local oscillator addition is shown in Figure 2. 
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Figure 2. Direction-finding System with Addition of Local Oscillator. 


This oscillator is used to downconvert the signals from the bridge-to-bridge 
frequency (between 156-162 MHz) to an intermediate frequency that is better suited for 
the filters in the system. As seen above in Section B, mixing the desired signal with its 
phase-shifted version enabled acquisition of the phase information. However, questions 
remained whether adding a local oscillator would still enable the phase information to be 
acquired. The following system of equations shows that the phase information is still 
able to be obtained, albeit with a more complicated system. For computational purposes, 
one assumes that the incoming signal has been adjusted to unit amplitude, and the local 
oscillator output has unit amplitude as well. The mixing of the incoming signal with the 


output of the local oscillator, yields output 


S, =cos(2nf.t)cos(2x(f, +z)t+8) 


(2.9) 
=0.5cos(22(2f, +z)t+0)+0.5cos(2zzt +6) 
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where the desired intermediate frequency is given by z and the phase difference between 
the incoming signal and the local oscillator is given by 8. Applying a low pass filter to 


remove the high-frequency (47 ft) components yields the following relationship: 
1 
LP{S,} = 5 008(2azt +8) (2.10) 


A similar method can be applied to the phase-shifted version of the desired signal. 
The following equations show the mixing of the incoming phase-shifted signal with the 


output of the local oscillator, yielding outputS,. As before, the desired intermediate 
frequency is given by z. 
S, =cos(2af,t+¢)cos(2x(f,+z)t+8) 


; : (2.11) 
S, = 5 cos| 20 (2, +2)t+9+0 |+—cos(2xz +0-9) 


Applying a low pass filter to remove the high-frequency (47 f.t) components 


yields the following relationship: 
1 
LP{S,} = —cos(2az +0 -4) (2.12) 


Thus from (2.12) one can see that the phase information is still preserved despite 
the addition of the local oscillator. Yet from (2.10) and (2.12), one can still see the local 


oscillator frequency remains in S,andS,. Thus, the signals need to be mixed again in 


order to achieve a function that solely depends on phase difference. As seen before, it is 
the phase information that is needed to obtain the Angle of Arrival. Using similar 


notation as before, a(t) is the mixed signal. 
a(t) =LP{S,}LP{S,} 
a(t) = +cos(2azt +8)x—cos(2nz1 +69) (2.13) 


a(t) = {005+ cos(4zt 20-9) 
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Applying a low pass filter to remove the double intermediate frequency 


(4zzt) components yields 
1 
LP {a(t)} = 0089 (2.14) 


This derivation shows that the solution is independent of the phase of the local oscillator. 
Like the result seen in (2.8), angle ambiguities arise when ¢ is used to calculate Angle of 
Arrival which will be detailed in III.B.1.c. Comparing (2.8) and (2.14), one can see that 
both methods yield results that are based solely on the phase difference. It was based on 
the result of (2.14) that design could proceed using the local oscillator method. Because 
a ship is a relatively slow-moving vessel, the phase difference will not be changing 
rapidly. Due to this fact, the output value of (2.14) is essentially a DC value. Thus, the 
phase can be applied to (2.4) to calculate the Angle of Arrival. 


Thus, one can see that is mathematically feasible to calculate the phase difference 
that occurs when one signal arrives at two different points. This information can then be 
used to calculate the Angle of Arrival from where the signal originated. Based on this 
foundation, work was undertaken to implement these findings in analog hardware, and 


this work is covered in the next chapter. 
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Hil. HARDWARE SOLUTION 


This chapter details the method of building the direction-finding system using 
analog hardware components. Simulation tools and techniques are described. Finally, 


the process of building the circuit boards is described, and results are presented. 


A. INTRODUCTION 


As shown in Chapter I, it proved mathematically feasible to determine the Angle 
of Arrival from the phase difference between two signals received at two separate 
antennas. Thus, an effective proof of concept was realized. The system was viable and 
capable of solving the problem investigated by this research. However, the equations and 
mathematical methods utilized were completely idealized. Critical assumptions that were 


made included: 
e Received signals were completely free from noise. 
e All mixing was achieved with zero loss. 
e All filters were completely ideal. 


The next step in the design process was to electronically simulate the entire 
system using established computer software before any hardware or components were 
purchased. This would ensure that the funds would not be wasted building a system that 
would not even perform in a computer environment. The simulation was performed on a 
system-wide level using Mathwork’s Simulink® [11] program, as well as on a system 


component level using Tonne Software’s Elsie™ [12]. 


B. SIMULATION TECHNIQUES 


This section describes the system-wide and lower level simulation techniques 


utilized in the hardware solution. 
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1. System-wide Simulation Using Simulink 


This section describes the system-wide simulation conducted without using a 
local oscillator. It then details the reason for the addition of a local oscillator to the 


system. 


a. System-wide Simulation without Local Oscillator 


oe 


From Mathwork’s corporate website, Simulink is described as “an 
environment for multidomain simulation and Model-Based Design for dynamic and 
embedded systems[11].” This software program was a logical next step in the design 
process, for a system-wide simulation was needed to test the proof of concept derived 
earlier. Simulink has a wide array of block libraries that lets the user design, simulate, 
and test both simple and complex systems. There are many user-customizable settings 
for the overall system, as well as for individual blocks within the system. These settings 
are able to be quickly changed, greatly enhancing the user’s ability to test the system. 
The system can be put through a wide range of testing conditions all from a personal 
computer. Simulink is also integrated with Mathwork’s flagship program, MATLAB®. 
Thus, the user is able to take advantage of MATLAB’s computational power as well as 
its programming capabilities. One also has the ability to run MATLAB scripts to develop 
companion algorithms to enhance the simulation, as well as to define signal, parameter, 
and test data [11]. It was these factors, along with Simulink’s prevalence in the 


engineering community, which made it a logical choice to model the desired system. 


Figure 3 shows the Simulink Model of the direction-finding system 


described in Chapter IL. 
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Figure 3. Initial Implementation of Direction-finding Model. 


The top Sine Wave block simulates the non-phase-shifted signal seen at 
one antenna. It has unit amplitude and a frequency of 156.8 MHz (bridge-to-bridge 
Radio Channel 16). White Gaussian noise is added to that sinusoid in order to simulate 
the noise added to a signal in a real-world environment. The combined signal is then 
passed through a Butterworth bandpass filter. The bandpass filter is needed so that the 
user can select which channel on the bridge-to-bridge frequency spectrum is of interest 
and to eliminate out of band signals, noise, and interference. The characteristics of the 
bridge-to-bridge frequency spectrum necessitate a very narrow bandpass filter. The 
bandwidth of each channel is only 50 kHz [1], which is much less than 1 percent of the 
channel frequency. This narrow passband and the close location of adjacent channels 
required a steep roll-off on either side of the passband. Thus, to only capture and 
analyze one channel, a high order bandpass filter is required. For the simulation seen in 


Figure 3, a Butterworth Filter with order 5 was used to analyze the system. 


When analyzing the above simulation, the passband was originally set 
strictly for Channel 16, with 156.8 MHz as the center frequency and 156.775 MHz and 
156.825 MHz as the 3 dB cut-off frequencies. This reflects the 50 kHz of bandwidth 
assigned to each bridge-to-bridge channel. With such a narrow passband, however, the 


system was not able to achieve an acceptable Angle of Arrival. The solution also did not 
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improve as the order of the filter was increased. Regardless of the filter order, the 
system was not able to determine the Angle of Arrival of the incoming signal. This was 
due to the fact that trying to achieve such a narrow passband resulted in unacceptable 


attenuation of the signal. 
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Figure 4. Bandpass Filter Centered at 156.8 MHz and 50 kHz Bandwidth. 


Figure 4 shows that with such tight parameters, there is almost 100 dB of 
attenuation at the center frequency. Due to this, the filter parameters had to be adjusted 


in order for the system to determine the Angle of Arrival. 


In order for the system to achieve an acceptable Angle of Arrival, the 
passband of the bandpass filter needed to be widened. The bandpass filter was then 
expanded to include the range of frequencies used by bridge-to-bridge radio, from 156.05 
MHz to 162.025 MHz [1]. The plot of the bandpass filter with these parameters is shown 


in Figure 5. 


16 


Transmission, dB 


i 





=5.833 


=11.67 





-17.5 


=23.33 





=29.. 17 



































150M 151.5M 153M 154.5M 156M 1S7.5H 159M 160.5M 162M 163.5M 165M 


Figure 5. Bandpass Filter Centered at 159 MHz and 6 MHz Bandwidth. 


By widening the passband, one can see that the passband attenuation has 
improved to approximately -5 dB. However, this has eliminated the ability of the filter to 
focus on one specific channel. If there is only one bridge-to-bridge channel being 
transmitted at a specific instant of time, then the system will be able to ascertain the 
Angle of Arrival of that signal. That scenario is unlikely, especially so in areas of high 
shipping density. In those situations, the more likely situation is that multiple ships will 
be having multiple conversations on a variety of bridge-to-bridge channels. When this 
happens, the system will not be able to track a specific signal, which is the overall 
purpose of the system. A bandpass filter, as shown in Figure 5, would be able to filter 
out other signals in the Very High Frequency band, such as KM radio broadcasts, air 
traffic control, and NOAA weather stations [1]. 


Because this research was intended to be used on United States naval 
vessels, which regularly operate in areas close to shore and in areas of high traffic, the 
above findings necessitated a re-design. A solution needed to be found that would enable 
the system to focus on one specific bridge-to-bridge channel. This would, in turn, 


eliminate all surrounding signals in the Radio Frequency spectrum. In addition to the 
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above requirements, the system would need to be able to be built and manufactured with 


off-the-shelf commercial components that are readily available and economical. 


b. System-wide Simulation Utilizing Local Oscillator 


As shown in Chapter II, Section C, it proved theoretically and 
mathematically possible to introduce a local oscillator to the system and still recover the 
desired phase information. The local oscillator would also address and solve the 
problems introduced in the previous section. Most importantly, utilizing a local oscillator 
in the system would eliminate the need for a bandpass filter. The system would become 
more complex with the addition of a signal generator and an additional stage of mixing. 
However, these additions would enable the system to be implemented using a 
combination of low-pass filters, which are significantly simpler to design, build, and 


implement. 


Utilizing a local oscillator would also enable the vessel operating the 
direction-finding system to focus on a specific bridge-to-bridge channel. Using the 
variables defined in Chapter II, one is able to design the system and first stage filters 
about a specific intermediate frequency z. Knowing the intermediate frequency, one can 


set the local oscillator frequency to( f+ z) . Thus, the system will only yield a coherent 


solution if the frequency of the incoming signal received at the antenna is equal to the 
local oscillator frequency, less the intermediate frequency. If they are equal, one is able 
to obtain the relationships defined in (2.10) and (2.12), shown again here for 


convenience. 


LP{S,} ==c0s (221 +6) (3.1) 


LP{s,} == cos(2rt+0-4) (3.2) 


Conversely, if the frequency of the incoming signal received at the 
antenna is not equal to the local oscillator frequency, less the intermediate frequency, the 


following relationship is seen: 
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S, =cos(27f,t)cos(2z(f, +z)t) 


S, = Scos[ 2a (f, te f,)t]+5c05[ 24 f, +z+f,)t| (3.3) 


where f.is the frequency of the desired bridge-to-bridge radio channel and fis the 
frequency of another signal, possibly another signal in the bridge-to-bridge frequency 
spectrum. When the incoming signal frequency matches the local oscillator frequency, 
less the intermediate frequency, a low pass filter is able to filter out high frequency 
components to achieve the relationships shown in (3.1) and (3.2). Applying a low pass 
filter will filter out the cosine and sine terms with a frequencies of ( fet f+ ae 
However, the cosine and sine terms with frequencies of ( tee Go f,) may not necessarily 
be eliminated. This is especially true if referring to two signals within the overall 
bridge-to-bridge frequency spectrum. Signals in that spectrum are separated by no more 
than 2 MHz, so the ( j.+2> f,) terms would not be eliminated by the low pass filter, 
with a 3 dB cut-off frequency. If those terms are not eliminated, one is not able to 
achieve relationships based solely on phase. One will never know f, due to the fact that 
it is an interfering signal, thus one will not be able to receive a coherent solution if the 
frequency of the signal at the antenna is not equal to the local oscillator frequency, less 
the intermediate frequency. One way to ensure that the ( Di ot oo f,) terms would be 
eliminated would be to put a tight bandpass filter on the front end of the system, with the 
passband set precisely at the desired f,. Another option would be to use a bandpass filter 
with 50 kHz bandwidth centered at the intermediate frequency. This filter would be used 
in lieu of the low-pass filters in the system. As stated earlier, such narrow bandpass 


filters proved infeasible to implement. These areas would benefit greatly from future 


research. 


The Simulink model shown in Figure 3 needed to be adapted to reflect the 
addition of the local oscillator to the system. The updated and final Simulink model of 


the direction-finding system is shown in Figure 6. 
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Figure 6. Simulink Model with Local Oscillator. 


As before, sinusoids represent the incoming bridge-to-bridge transmission. 
For the simulation, they were given unit amplitude and set at 156.8 MHz, representing 
Channel 16 on bridge-to-bridge radio. In this final model, a third Sine Wave block has 
been added to represent the local oscillator. The oscillator frequency was set to 167.5 
MHz (156.8 + 10.7 MHz). The intermediate frequency (IF) was chosen to be 10.7 MHz, 
due to the prevalence of this frequency as an IF in FM radio circuits [13]. Because this 
system is meant to be an economical solution to the direction-finding problem, using 
technical specifications that are common in the commercial world will lead to lower 


implementation and fielding costs of the system. 


Additive white Gaussian noise is added to the incoming signal and its 
phase-shifted version in order to simulate the noise encountered in a real-world 
environment. These signals are then multiplied by the local oscillator output and passed 
through a low-pass filter. In Simulink, these filters were modeled by third-order 
Butterworth low-pass filters. The cut-off frequency was set to the IF frequency of 10.7 
MHz, which will pass the information contained at that frequency and below but 
eliminate the information at the original carrier and higher frequencies. After this first 
round of filtering, the system has the relationships seen in (2.10) and (2.12). These 


signals, S, and S,,, are then mixed together and passed through another low-pass filter. 
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This low pass filter is designed to filter out the intermediate frequency and leave an 
essentially DC signal that is a function of the phase difference. This low-pass filter was 
given a cut-off frequency of 1 kHz. This value was chosen as it is significantly lower 
than 10.7 MHz, but it also accounts for the fact that in a real-world environment the final 
signal will not be purely DC. After filtering, the signal is multiplied by a constant value 


of eight. After that multiplication, the signal is now equal tocos¢ . 


c. Requirement for Two Systems 


After the signal passes through the constant gain of eight, calculations 
must be performed to transform the calculated phase into an Angle of Arrival. The first 


step in this process is taking the inverse cosine of the signal, which should yield the phase 
difference. However, the inverse cosine is only defined on the interval [0, z] . What this 
means is that for every value of cos¢ , there are two possible phases on [0,27]. Thus, the 
system must then be split into two to take into account the two possible phases. 


The system for translating the phase difference into Angle of Arrival is 


shown in Figure 7. 
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As seen above the system is split into two arms to account for the two 
possible solutions to the inverse cosine operation. The signal in the bottom arm is 


multiplied by negative one to account for angles between|z, 2r]. As shown previously 


in (2.4), the relationship between phase difference and AOA is: 


AOA = arcsin Cuma (3.4) 
360-B:- f. 


Converting (3.4) into radians gives the following relationship: 





408-=asin{ eg (3.5) 
22 f.B 


When designing the Simulink model, it was assumed that the baseline 
distance between the two receiving antennas was equal to (44). where / is the 
wavelength of the desired bridge-to-bridge signal. For bridge-to-bridge radio channel 16, 
yy, is approximately 0.96 meters. As stated in Ch. II, accuracy will improve as the 


length of the baseline, relative to the transmitted signal’s wavelength, is increased. This 


is true up to a baseline equal to one-half wavelength. After that point, the phase 


difference changes more than 360° as the AOA moves from +90°to —90°[9]. Thus, 


assuming b= (44), (3.5) simplifies to: 


AOA = arcsin (4) (3.6) 
1 


After performing the inverse cosine operation, one has two possible 


solutions for ¢, one solution between [0,7] and the other between [7,27]. One must 


multiply both values of ¢ by (Y) to achieve the quantity (%) . Because 


sin( AOA) = (4) (3.7) 
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the final step to obtaining the Angle of Arrival is to take the inverse sine of the signals in 


both arms. However, like the inverse cosine, the inverse sine is only defined for a limited 


range of angles. The inverse sine is defined over the range |-2.2] . Because of this 


property, the two arms must then be split into four branches to take into account the range 


of angles EE This will account for all possible Angles of Arrival. The 


mathematical operations to achieve all four possible Angles of Arrival are detailed below. 
For the first branch, after the arcsine operation is performed, the signal is multiplied by 


(180/) to convert the signal from radians to degrees. By performing this operation, the 
system is able to calculate Angles of Arrival between [ 0°. 90" | . The relationship 
between the output of the first branch and the third branch 1s: 


Third Branch Output = (180 — First Branch Output ) (3.8) 


This operation allows the system to calculate Angles of Arrival between [:90°,180° }. 
The second arm (in which the output of the inverse cosine was multiplied by negative 
one) is multiplied by (1807) to convert the signal to degrees. That signal is then added 


to 360 in order to obtain a positive Angle of Arrival; this produces the output of the 


second branch. By this operation, the system is able to calculate Angles of Arrival 


between [ 270° : 360° | . Finally, the relationship between the second and fourth branch is: 
Fourth Branch Output = (180 - Second Branch Output). (3.9) 


This operation takes into account Angles of Arrival between [ 180°, 270" ]. While these 


operations add complexity to the system, they are necessary to ensure the system is able 


to calculate Angles of Arrival from the entire range of possibilities [ 0°.360° | : 
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The purpose of this research was to develop a system that could be utilized 
onboard United States vessels with minimal user input. It should be a system from which 
the user can receive a clear solution and then act accordingly. If the system gives four 
possible answers, the user is left to decipher which answer is correct. In an environment 
in which traffic is sparse, this would not be a difficult task. One could simply look in the 
four directions and see clearly which of the four angles had other vessels on axis. 
Process of elimination would lead the watchstander to know which ship was transmitting. 
However, if traffic is heavy, there could be ships at each of the four possible Angles of 
Arrival. The watchstander would have absolutely no way of knowing which ship was 
transmitting. The system would need to be changed in order to eliminate the ambiguity a 


watchstander would encounter. 


The solution that was determined was to introduce another pair of 
antennas offset from the original pair at a pre-determined angle. Thus, each pair of 
antennas would produce four solutions, and the solutions from Antenna Set 1 would be 
compared to the set of solutions from Antenna Set 2. By comparing solutions, one could 
find commonalities and determine the true Angle of Arrival of the transmitted signal. 
Different offset angles were experimented with to determine if a single, coherent solution 
could be obtained. The first offset attempted was a 90-degree offset of the second pair of 
antennas. Unfortunately, this amount of offset was unable to narrow down the solution 
set and four possible solutions remained. After multiple attempts with a variety of angle 
offsets, the best solution was obtained with an offset of 45 degrees. While this solution 
was not able to determine one single solution, it did narrow the solution set down from 
four possible solutions to two possible solutions. Each of the two antenna systems was 
analyzed at every Angle of Arrival between 0 and 360 degrees. The solution sets were 
then compared at every angle and certain patterns were noticed. At certain angles, for 
example, Solution Set 1 from antenna system | matched Solution Set 1 from antenna 
system 2. At the same angles, Solution Set 4 from antenna system 1 matched Solution 


Set 4 from antenna system 2. What made this finding important was that these solution 
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sets only matched at certain Angles of Arrival. Thus, by using a system of logic 
statements one was able to narrow the solution set down from four to two. Table 1 shows 


the input to the logic statements used to narrow the solution set. 


True Angle of Arrival (degrees) |System 1 |System 2 |System 1 
AOA AOA4 
AOA3 AOA4 



































Table 1. | Method for Determining Solution Set. 


For example, if Antenna System | returns the solutions in Table 2: 























AOA1 30 
AOA2 330 
AOA3 150 
AOA4 210 





Table 2. | Example Solution Set from Antenna System 1. 


and Antenna System 2 returns the solutions in Table 3 (after subtracting out the angle 


offset of the second antenna system), 














AOA1 30 
AOA2 240 
AOA3 60 
AOA4 210 














Table 3. | Example Solution Set from Antenna System 2. 


one can clearly see that AOAI from Antenna System 1 matches AOA1 from Antenna 
System 2, and AOA4 from Antenna System 1 matches AOA4 from Antenna System 2. 
Additionally, no values from either Antenna System’s AOA2 or AOA3 are the same, so 
the solution set is down to two possible values: 30 degrees or 210 degrees. While this 


solution is not ideal for the watchstander, it has given him or her only two possible 
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solutions, separated by 180 degrees. The inability of this system to yield one definitive 
solution is a weakness in the system. Further research is recommended in this area to 


ascertain methods or algorithms to solve this problem. 


Ze Lower-level Simulation Using Elsie 


Simulink effectively confirmed that the mathematical derivations could be 
implemented on a system-wide level and generate a set of two possible solutions. 
However, Simulink is still a highly idealized system. The Product blocks in Simulink 
were replaced by analog mixers when the system was actually built. The multiplication 
in Simulink is ideal, but analog mixing results in some power loss. Likewise, a perfectly 
ideal filter does not exist in the real world. When building an analog filter, achieving the 
exact cut-off frequency one desires is an unrealistic expectation. The next step in the 
design process was performing a more realistic simulation of the filters contained within 
the system. The software program chosen to perform this filter simulation was Tonne 


Software’s Elsie, an electrical filter design and analysis program. 


From Tonne Software’s website, Elsie is described as a “32-bit Windows 
electrical filter design software nicely written to help engineers design and analyze 
lumped-element filters in the audio through microwave range.” It allows the user to 
design lowpass, highpass, bandpass, and bandstop filters in many families of filters, 
including Butterworth, Chebyshev, Cauer, Bessel, and Gaussian. The student version of 
Elsie, which was utilized during this research, allows the design of filters up to seventh 
order. The program contains significant plotting capability, allowing thorough analysis 
of the designed filter. Once the ideal filter is designed, the program allows the user to 
substitute real-world component values into the system to see how the real-world filter 
differs from the theoretical version. Elsie proved a very capable program for linking the 


theoretical design of the filter to the building of physical filters [12]. 


The first step in using Elsie was to design the first stage of low-pass filters. As 
noted in the previous section, this filter was to pass frequencies up to the intermediate 
frequency of 10.7 MHz. The first attempt at the first-stage filter was that of a third-order 


Butterworth low-pass filter. The frequency response is shown below in Figure 8. 


26 


Transmission, dB 


0 












































—_—_—. | 
@ 

-5.833 
-11.67 
-17.5 
-23.33 
-29.17 

=35 

1K 2.0009M 4.0008M 6.0007M 8.0006M 10.001M 12M 14M 16M 18M 20M 
Figure 8. Ideal First Stage Low-Pass Filter. 


The marker shown on Figure 8 is located at 10.7 MHz, with a transmission loss of 
3 dB. This frequency response can be obtained with the third-order Butterworth filter 


shown in Figure 9: 
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Figure 9. Ideal First Stage Low Pass Filter Schematic. 
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As one can clearly see in Figure 9, the capacitor and inductor values are purely 
theoretical. It is not possible to purchase such components accurate to the third and 
fourth decimal point commercially. Elsie was then used to convert these component 
values into components that could be purchased commercially. Figure 10 shows the 


schematic generated by Elsie after conversion to the commercial components. 
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Figure 10. First-Stage Filter with Commercially Available Components. 


Changing the component values naturally changed the frequency response of the 


filter. The changed frequency response is shown in Figure 11. 
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Figure 11. Frequency Response of First Stage Filter with Commercially Available 
Components. 
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By changing the component values, the transmission loss has increased from 3 dB 
to 3.17 dB at 10.7 MHz. This corresponds to a loss of 52 percent of the power after 
passing through just the first stage filter. Alternate filter designs were then explored with 
the possibility of lessening the power loss after the first-stage filter. After multiple 
iterations, the final design was reached with a cut-off frequency of 20 MHz. While this is 
greater than the intermediate frequency of 10.7 MHz, it is still significantly lower than 
the bridge-to-bridge frequency range, which begins at 156.05 MHz [1]. Staying 
significantly lower than the bridge-to-bridge frequency range ensures that the filter is 
capable of eliminating the information contained in the higher frequencies, as shown in 
(2.10) and (2.12). The frequency response of the final design of the first-stage filters is 


shown in Figure 12. 
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Figure 12. Frequency Response of Final Design of First-Stage Filters. 


As in the previous figures, the marker is set at 10.7 MHz. This design led to 
significant improvement in the power loss, as there was just 0.91 dB of power loss at 10.7 
MHz. This is equivalent to just 19 percent power loss, while still keeping the filter 
simple for construction purposes. The above performance was achieved with a third- 


order filter. Figure 13 shows the final filter schematic of the first-stage low-pass filters. 
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Figure 13. Circuit Schematic of Final Design of First-Stage Low-Pass Filters. 


Finding a solution that could be accomplished with a third-order filter was 
important. Not only would the system be easier to implement, it would also be more 
cost-effective. Additional complexity would not only make building the system more 


difficult, but adding additional components would significantly escalate cost. 


After the first-stage filters were designed, the next step was to design the second- 
stage filter. As noted in the previous section, this filter was to pass frequencies from DC 
up through | kHz. The first attempt at the second-stage filter was that of a third-order 


Butterworth low-pass filter. The frequency response is shown below in Figure 14. 
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Figure 14. Ideal Second-Stage Low-Pass Filter. 


The marker shown on Figure 14 is located at 1 kHz, with a transmission loss of 
3.06 dB. The following filter design can be accomplished with a third-order Butterworth 


filter using the following schematic: 
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Figure 15. Ideal Second Stage Low Pass Filter Schematic. 


31 


As one can clearly see in Figure 15, the capacitor and inductor values are purely 
theoretical. It is not possible to purchase such components accurate to the fourth decimal 
point commercially. Elsie was then used to convert these component values into 
components that could be purchased commercially. Figure 16 shows the schematic 


generated by Elsie after conversion to the commercial components. 
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Figure 16. Second-Stage Filter with Commercially Available Components. 


Changing the component values naturally changed the frequency response of the 


filter. The changed frequency response is shown in Figure 17. 
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Figure 17. Frequency Response of Second Stage Filter with Commercial Parts. 


By changing the component values, the transmission loss has increased from 3 dB 
to 3.27 dB. This corresponds to a loss of 53 percent of the power after passing through 
the second stage filter. Alternate filter designs were then explored with the possibility of 
lessening the power loss. After multiple iterations, the final design was reached with a 
cut-off frequency of 1 kHz. The frequency response of the final design of the first-stage 
filters is shown in Figure 18. The rise in the frequency response at approximately 6 kHz 
is due to the fact that a desired inductor was not commercially available and near- 


substitute component had to be used in the circuit construction. 
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Figure 18. Frequency Response of Final Design of Second-Stage Filter. 


Figure 19 shows the final filter schematic of the second-stage low-pass filter. 
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Figure 19. Circuit Schematic of Final Design of Second-Stage Low-Pass Filter. 
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As before, finding a solution that could be accomplished with a third-order filter 
was important. Not only would the system be easier to implement, it would also be more 
cost-effective. Additional complexity would not only make building the system more 


difficult, but adding additional components would significantly escalate cost. 


Thus, the final system could be built utilizing only third order filters and 
commercial off-the-shelf components. By accomplishing these objectives, the next step 


of building the system would be even more manageable. 


C. BUILDING THE HARDWARE SOLUTION 


The next step in the design process was to build the circuits that would be 
required to implement the direction-finding system. Research was conducted into 
various electronics companies to ensure that the component values determined in the 
Elsie simulations were available for purchase. The companies, on the whole, carried the 
same components and in the same magnitudes. However, they differed in the availability 
of components, the lead time required for purchase, purchase minimums, and other 
factors. A major decision was whether to use through-hole components or surface- 
mounted components. Another consideration was the compatibility of the components 
with the expected test equipment. The capacitor values used in the first-stage filters are 
in the picofarad range, and the capacitors used in the second-stage filter are in the 
microfarad range. With such small values, the issue of stray capacitances was very 
important. Ultimately, surface-mounted components were chosen based on the fact that 
there would be much less stray capacitance than if through-hole components were used 
[14]. The connectors on the edge of each circuit board needed to be SMA (SubMiniature 
version A)-compatible, so they would be compatible with the test equipment. This 
requirement limited the number of connectors available for purchase. Ultimately, the 
most cost-effective connectors that fulfilled the design requirements were through-hole 
components. The capacitors, inductors, and resistors were able to be purchased as 


surface-mounted components. 
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Before the components could be purchased, the circuit boards on which they 
would be mounted needed to be designed and purchased. The company PCB123™, 
whose parent company is Sunstone Circuits®, was selected to manufacture these circuit 
boards [15]. PCB123 was ideal for a number of reasons. The design software for the 
circuit boards is available as a free internet download, with free customer and technical 
support. Once one is familiar with the software, the user is able to design the circuit 
boards with relative ease. The software comes with a comprehensive part library. This 
was advantageous because the boards could be designed with the exact components that 
were determined with the Elsie software. While the circuit boards required for this 
research needed only two electrical layers, PCB123 allows up to six layers. Another 
important feature was that after design was complete, the designs were sent electronically 
to PCB123 for manufacture. Not only would PCB123 manufacture the boards, the 
company would perform electrical checks afterwards to determine that the boards were 
manufactured correctly [15]. In the design process, the user is able to specify the 
locations and layout of all components, the size of the components, and the electrical path 
connecting the filters. Figure 20 shows the PCB schematic that incorporates both first- 


stage low-pass filters and the second-stage low-pass filter. 
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Figure 20. PCB 123 Circuit Schematic. 
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Figure 20 shows the theoretical schematic of the direction-finding system on the 
circuit board. The PCB123 software converts the schematic to a “layout” model, where 
one can manipulate the physical location of components and electrical paths. 


1st Stage Low Pass 2nd Stage Filter 
Filter | 





Figure 21. PCB 123 Layout Model. 


Figure 21 shows the layout model generated by PCB123. The boards chosen 
were three inch by five inch circuit boards. There were six SMA connectors per board, 


which are the large square objects seen near the borders in Figure 21. 


The direction-finding system operates as follows: The incoming signal is 
multiplied by the local oscillator using an analog mixer, with the output entering the top 
first-stage low-pass filter from the left. The incoming phase-shifted signal is multiplied 
by the local oscillator using an analog mixer, with its output entering the bottom first- 
stage low-pass filter from the left. Both signals pass through their respective low-pass 
filters and then exit the board through the SMA connectors located in the middle of the 
circuit board. Those two signals are then multiplied together, by way of an analog mixer, 
and proceed to the second-stage low-pass filter. The output of the mixer re-enters the 
board through the top of the second-stage filter. The signal then passes through the 
second-stage low-pass filter and the DC voltage is read by a voltmeter for analysis. 
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D. HARDWARE SOLUTION RESULTS 


After the circuit boards were received, the electrical components were ordered. 
The components were then soldered onto the circuit boards. Each antenna system 
requires one circuit board; thus, a minimum of two boards are required for the direction- 
finding system. Four were created, however, to increase the probability of having a 
functioning circuit board. After construction, the circuits were tested to ensure the filters 
met the design criteria. All four circuit boards were independently tested, and all four 
circuit boards yielded the same testing results. As stated in previous sections, the first- 
stage filters were built with a cut-off frequency of 20 MHz. Thus, there should be 3 dB 
of attenuation at 20 MHz. When tested, however, the first-stage filters experienced 3 dB 
of loss at 15 MHz. This was not a problem, though, for the desired intermediate 
frequency of 10.7 MHz was still well within the passband. In essence, the filter was 


tighter than designed but would have no negative effect on overall system performance. 


Like the first-stage filters, the second-stage filter on each circuit board was tested 
independently. Each of the second-stage filters on the four circuit boards all had nearly 
identical performance. The second-stage filters were designed to have 3 dB of 
attenuation at | kHz. Upon testing, it was determined that there was actually 4 dB of 
attenuation at 1 kHz. While this is worse performance than originally designed, the 
signal of interest is basically DC, so the decreased performance at | kHz did not have an 


overall negative effect on system performance. 


After verifying each circuit board’s individual performance, the overall system 
was constructed. A signal generator was set at 156.8 MHz to simulate bridge-to-bridge 
radio Channel 16. The signal was also frequency modulated, to more realistically 
simulate an actual bridge-to-bridge radio signal. Another signal generator was used to 
simulate the local oscillator; this generator was set to 167.5 MHz, which reflects the 
bridge-to-bridge frequency plus the intermediate frequency. The system was then set-up 
as described in the previous section. The signal from the first signal generator was split, 
with one output unchanged and the second output passing through a Sage Laboratories® 
Model 6708 phase shifter [16]. This phase shift represented the phase shift that would be 


caused by the bridge-to-bridge signal arriving at the second antenna. Both signals were 
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multiplied by the local oscillator signal with a Mini Circuits® ZLW-1 analog mixer [17]. 
The local oscillator used was a Hewlett Packard 8341B Synthesizer Sweeper [18]. Both 
signals were then passed through the first-stage filters on one circuit board. After the 
signals exited the first-stage filters, they were mixed together with an analog mixer and 
then passed through the second-stage filter. The output of the second stage filter was 
then sent to a voltmeter for analysis. A picture of the overall testing system is shown in 


Figure 22. 





Figure 22. Hardware Solution Testing System. 


The phase shifter used in the experiment varied the phase from zero to 180 
degrees. The phase was changed by use of a knob on the phase shifter. As the knob was 
turned, a change in phase was introduced. The range of the phase shifter was used to 
determine how many degrees of phase shift were introduced by each turn of the knob. 


After the calibration of the phase shifter was verified, the experiment was then conducted 
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multiple times by comparing the voltage output of the system to the phase difference 
between the two original signals. Interestingly, throughout multiple trials, a clear linear 
relationship was observed. This was unlike what was seen in (2.14) and shown again 


below in (3.10) 
LP {a(t)} = cos (3.10) 


where a trigonometric relationship was derived. This was an unexpected benefit, as the 
periodic nature of trigonometric functions can yield solutions that are possible from more 
than one phase difference. The linear solution yields one unique solution per phase 
difference. Further research should be conducted to determine why the linear 
relationship was experienced, rather than the expected trigonometric relationship. The 
relationship is seen in the figure below, with the experimental data marked with circles 


and a best-fit line superimposed. 
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Figure 23. Relationship Between Phase Difference and Voltage Output. 
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One can clearly see the linear relationship between the phase difference and the 
voltage output. Using the equation for the best-fit line, the relationship between phase 


and voltage output is: 
deg 
phase(deg) =| —62.5 ray I voltage(mV )—615.625(deg ) (3.11) 
m 


This phase information is then entered into a MATLAB program that uses the 
equations derived in Chapter II to arrive at a solution set of possible Angle of Arrivals for 
Antenna System 1. Antenna System 2 is set up in the exact same way as Antenna System 
1, but with a pre-determined angle offset. Comparing the solutions obtained from 
Antenna System 1 with the solutions obtained from Antenna System 2, one is able to 


determine the Angle of Arrival of the incoming signal. 


Thus, the hardware solution has solved the problem for which the research was 
conducted. It has proved that a functioning system can be built with analog components 
that can accurately determine the Angle of Arrival of an incoming bridge-to-bridge signal 
if the phase difference between a single signal arriving at different locations is known. 
As the hardware solution was being explored, a software solution was being explored 
concurrently. Obtaining a solution in software would allow two different solutions to the 
direction-finding problem, which could be compared and analyzed to determine the best 


solution. 
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IV. SOFTWARE SOLUTION 


This chapter details the method of building the direction-finding system using 
software-defined radio. The software projects utilized are introduced and described, as 
well as the components used in the system. Finally, the process of setting up the 


experiment is described, and results are presented. 


A. INTRODUCTION 


As seen in Chapter II, it proved mathematically feasible to determine Angle of 
Arrival from the phase difference between two signals received at two separate antennas. 
The findings in Chapter III verified that a solution was possible using hardware analog 
components. Research was conducted concurrently to determine if a direction-finding 


solution could be obtained using software. 


Rapid improvements in technology, have led to many functions traditionally 
found in hardware to now be implemented in software. Software-based solutions are 
often preferable as they can be cheaper, more flexible, and achieve higher performance. 


It is these advantages that led to the development of the software radio [19]. 


Ideally, the software-defined radio receiver would have an analog-to-digital 
converter (ADC) attached directly to an antenna. The ADC would digitize the analog 
signal, and then computer software would process all the results and perform any needed 
applications or calculations. However, due to the desired frequencies of the signals of 
interest, oftentimes a downconverter is required to bring the signal to a lower 
intermediate frequency or baseband before digital sampling can occur [10]. This device 
utilized in this thesis to do both downconversion and analog to digital conversion is the 
Universal Software Radio Peripheral (USRP) [20]. The first step in developing a digital 
solution to this problem was finding a software program that would be able to 
communicate with the USRP as well as be capable of implementing the direction-finding 


system. 
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B. GNU RADIO SOFTWARE APPLICATION 


Two of the most popular open-source software projects for software-defined radio 
currently used are the GNU radio project and the Open Source SCA Implementation 
Embedded (OSSIE) project. The OSSIE project is a Virginia Tech University initiative 
to provide a Software-Defined Radio platform for the development of waveforms based 
on the Software Communications Architecture specifications under the Joint Tactical 
Radio System program. OSSIE, originally released in 2004, was developed primarily for 
research and education in software-defined radios. Additionally, the Naval Postgraduate 
School has worked with Virginia Tech in development of laboratory exercises. The 
OSSIE project is written in C++, and the focus on development is on the Linux operating 
system. However, there has been initial research performed on implementing OSSIE on 


other operating systems [21]. 


The other major Software-Defined Radio Project, and the project used in this 
research, is the GNU Radio Project. Like OSSIE, GNU Radio is a free, open-source 
Software-Defined Radio development platform. It contains low-level programs and 
different processing blocks to implement these radios. GNU Radio contains many 
libraries for modulation, error-correcting codes, and signal processing. GNU radio code 
is primarily written using the Python programming language, while the supplied, 
performance-critical signal-processing path is implemented in C++ using processor 
floating point extensions. The strength of GNU radio is that the waveforms transmitted 
and received by the radio are defined by software. This gives GNU radio and other 
software-defined radio systems their biggest advantage—reconfigurability. Multiple 
expensive radios no longer need to be purchased; a single generic radio can be purchased 
and then be programmed with signal processing software to perform the desired function. 
The same generic radio can then be reprogrammed if the desired functionality and 


purpose of the radio system changes [3]. 


For this research, the digital solution to the direction-finding problem was 
explored using GNU radio software. However, only half of the functionality of the 
software was utilized. This research did not require the use of the transmission capability 


of the radio. No waveforms needed to be developed, nor transmitted. The solution being 
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researched was a passive system, which required only the receive capability of the GNU 
radio software. The bridge-to-bridge radio signal of interest is an FM-modulated analog 
signal. Because of this property, an analog-to-digital converter (ADC) was needed that 
could convert the VHF signal to a digital signal that can be analyzed and processed by the 
GNU radio software. The Universal Software Radio Peripheral contained the ADC 


utilized in this research. 


C. UNIVERSAL SOFTWARE RADIO PERIPHERAL 


The Universal Software Radio Peripheral (USRP), developed by Matt Ettus, is a 
low-cost software radio device that connects to a personal computer using a USB 2.0 
connection. It consists of a motherboard, which performs the analog-to-digital 
conversion of the baseband signal, and up to four daughterboards, which convert the 
original analog signals to the desired Intermediate Frequency (IF). The motherboard 
contains a Field Programmable Gate Array (FPGA) for high speed signal processing. 
The motherboard is shown in Figure 24 [20]. 
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Figure 24. USRP Motherboard. (From [20]). 


Each daughterboard can be used for a different frequency band, making the USRP 
an extremely flexible and useful device. Each motherboard can support up to two 
transmit daughterboards and up to two receive daughterboards. The basic receive 


daughterboard is shown in Figure 25 [20]. 
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Figure 25. Basic Receive Daughterboard. (From [20]). 


Utilizing the various daughterboards, the USRP has an overall range of DC to 5.9 
GHz, making it a very useful tool for the radio frequency spectrum [20]. The USRP 
works with GNU radio, an open-source framework for the creation of software-defined 


radios [3]. 


1. USRP Receive Capabilities 


The USRP daughterboard is responsible for transmitting the separated I and Q 
channels of the incoming complex signal to the motherboard and converting the received 


real signal to a complex signal at an intermediate frequency. The USRP contains four 


high-speed 12-bit analog-to-digital converters with sampling rates of up to 64x 10° 
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samples/second [20]. This gives a theoretical upper frequency limit of 32 MHz for the 
intermediate frequency before aliasing artifacts are introduced [22]. Two ADCs are 
required to convert the complex analog signal to a digital signal. One ADC is required 
for the I channel, and one ADC is required for the Q channel. The full voltage range of 
the analog-to-digital converters is 2 volts peak-to-peak, and the input is 50 ohms 
differential. This translates to 40 mW, or 16 dBm [22]. In the case of a weak received 
signal, the daughterboards can utilize programmable gain amplifiers before the ADCs to 
amplify the input signal up to 20 dB. The USRP also contains four digital-to-analog 
converters for transmission purposes, but as noted before, the transmission capabilities of 


the USRP were not required for this research [20]. 


After sampling, the motherboard FPGA processes the samples, essentially 
performing high bandwidth mathematical operations to reduce the data rate to an amount 
that is compatible with USB2.0 technology [22]. The standard FPGA configuration 
consists of digital down converters combined with low-pass filters. The decimation 
process is accomplished by low-pass filtering followed by downsampling. The USRP 


performs the filtering by means of cascaded integrator-comb filters [22]. 


Aliasing is a condition that can possibly result from digital sampling. It occurs 
because frequencies that differ by a multiple of the sampling rate are indistinguishable 
from each other after sampling [23]. The relationship between the angular analog 


frequency ©), and the angular digital frequency «@, is 


Q, a 2nf, 
i (4.1) 


Ss 


@ =Ogt, = 27 


In the above equations, f. is the frequency of the analog signal, and f, is the sampling 
frequency. One can clearly see that the digital frequency is dependent on the sampling 
frequency. When the FPGA performs the decimation, it is effectively lowering the 
sampling frequency, which changes the digital frequency. Within the digital frequency 


domain, any frequency @,and @,are indistinguishable from each other given 


that @, =@,+k2z, where kis any integer [10]. Frequencies above the Nyquist rate are 
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indistinguishable from the frequencies below it, which is the basis for aliasing. Because 
the Nyquist frequency, which is defined as one-half of the sampling rate, corresponds to a 


digital frequency ofz, a digital signal’s bandwidth must be entirely within the range 


[-2, x] to prevent aliasing [23]. Using (4.1), it follows that f. must be no greater than 


4 to ensure aliasing will not occur. 


2. USRP Transfer Capabilities 


The theoretical maximum data transfer rate of a USB 2.0 connection is 480 
Megabit/second [24]. In practice, however, this theoretical rate cannot be achieved; a 
more practical rate to be considered is approximately 256 Megabit/sec or 32 
Megabytes/sec [10]. Because each complex sample that is transferred from the USRP to 
the personal computer is 16 bits per channel (I and Q), the total size of a sample is 32 


bits. Thus, the practical sample rate utilized in this research is: 


« bits : lsample 





R, = 256 x 10 ; 

sec 32 bits (4.2) 
R =8x 10° samples 

: sec 


Utilizing a lower sample rate can prove more beneficial for computer processing as it will 


require less storage. This will be detailed in the next section. 


D. SOFTWARE SOLUTION EXPERIMENTATION 


After the capabilities of GNU radio and the associated USRP were known, the 
next step in the research process was to integrate them through programming to 
demonstrate the theory developed in Chapter II. Unlike in the analog solution, only one 
signal generator was required for this experiment. This is due to the fact that the single 
signal generator can be used to simulate the bridge-to-bridge signal. A local oscillator is 
not used in this solution, eliminating the need for the second signal generator utilized in 
the analog solution. To simulate the bridge-to-bridge signal arriving at the two shipboard 
antennas, a phase shifter was utilized to implement a phase shift into the transmission. 


To simplify the experiment and to provide consistency between the hardware and 
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software solutions, the same Sage Laboratories® Model 6708 phase shifter [16] was used 
in the software solution as was used in the hardware solution. By using the same phase 
shifter, equivalent test conditions could be established for both the hardware and software 
solutions. Thus, using the same components in both experiments will allow for more 
valid comparisons with the results. A block diagram of the software solution experiment 


is shown in Figure 26. 
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Figure 26. Software Solution Experiment Block Diagram. 


For the experiment conducted with the software-defined radio, the output of the 
signal generator and its phase-shifted version were passed through an analog mixer 
before entering the receive daughterboard of the USRP. Ideally, the mixing would be 
performed in software along with all signal processing. Time constraints prevented the 
implementation of the mixer in software, but that is recommended for further research on 
this topic. Once, the output of the mixer entered the USRP, a program was needed to 
acquire the data and convert it to a format that could be analyzed and yield constructive 


results. The signal entering the USRP should ideally first be passed through a bandpass 
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filter, which would enable selection of a specific channel in the bridge-to-bridge 
frequency spectrum. Time constraints prevented implementation of this filter, but further 


work into this filter implementation is highly recommended. 


As stated earlier, the GNU Radio Project contains numerous libraries for 
modulation, error-correcting codes, and signal processing. Along with the various 
blocks, the software download includes sample programs that can be used or modified to 
suit the user’s purpose. That is a major advantage to open-source software-defined radio. 
By letting many users use, change, and improve the software, the project has the ability to 
grow at a fast rate and more economically than proprietary competitors. The sample 
program that I used as the base for my signal capturing program is titled usrp_rx_cfile.py. 
This program, written in the python programming language, reads samples from the 
USRP and then saves the data into a user-specified file format. The user is able to 


specify many different factors into how the program will capture the data from the USRP. 


The first option lets the user select from which of the two daughterboards to 
collect data. As stated earlier, the USRP is able to contain up to two received 
daughterboards and up to two transmit daughterboards. There was only one 
daughterboard on the USRP used for this research, so this program specification was 
clear. This research utilized the BasicRX daughterboard [25] in the motherboard’s 
“RXA” position. However, situations could arise where both daughterboards are 
receiving data, but the user only wishes to analyze the data from one daughterboard, but 
not the other. One is also able to set a limit on the number of samples to collect, any 
desired gain, as well as an option to produce 8-bit samples instead of the traditional 16-bit 
samples. The command line specification that had the greatest impact on this research 
and required the most experimentation was the desired decimation. One would normally 
desire that the decimation rate be as small as possible to obtain as many data samples as 
possible during a given time period. However, in the shipboard environment envisioned 
by this research, the relative phase change between vessels changes slowly. This will 


allow the user the ability to use a higher decimation rate. As shown in (4.2), the 


. E samples 
maximum sample rate for the USB2.0 connection was8x 106 SAP** Because the 
sec 
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0° samples 


analog-to-digital converters in the USRP operate at 64x 1 , this corresponds to 





sec 
a decimation rate of 
Re Rosep 
dec Rea 
P (4.3) 
_ 64x10? | 
“8 x10° 


Thus, eight is the minimum decimation rate possible for the USB2.0 connection. 
If the decimation rate were less than eight, the data rate would be too high for the USB2.0 
to transfer to the personal computer [24]. Additionally, storage requirements of the 
personal computer played in role in determining the proper decimation rate for the 
experiment. Because each complex sample contains 32 bits of data (16 each for the I and 
Q channels), the storage requirement for just one second of data collection is shown 


below [22]: 





i 1 
Storage -( ozeits |[: x 19° Supers =| 
sample sec 


Storage = 256 x 10° eS 2405 10° bytes (4.4) 
sec sec 


MB 
Storage = 32 —— 
sec 


This is a substantial amount, especially considering that the personal computer 
has purposes other than just storing data. It must have enough memory for the signal 
processing software, as well as its other operating requirements. In a high traffic 
shipping environment, the shipboard user may want to use the program numerous times 
in a short amount of time. Even if the ship is only tracking one vessel, it will want to use 
the direction-finding system on a periodic basis, regularly updating its position. Thus, the 


storage requirement was a major design specification. 


a2 


The chosen decimated rate was also affected by the USB connection. As shown 
in (4.2), the maximum sample rate for the USB2.0 connection was calculated to 


0° samples 


be8x 1 . USB2.0 technology has incorporated buffers into its protocol. These 


sec 

buffers are intended to “compensate for a difference in data rates or time of occurrence of 
events, when transmitting data from one device to another [24].” However, if the 
decimation rate was too low, the system experienced overflows of the USB buffer. 
Experimentation was conducted to determine the effect of different decimation rates. A 
system trade-off had to be conducted in order to choose between a low decimation rate 
giving a greater number of samples versus a high decimation rate which was preferable 
for storage and USB buffer requirements. Ultimately, a decimation rate of 16 was 
chosen, for it was the lowest decimation rate which yielded coherent results. This 


lowered the USB sampling rate to: 


_ 64x10° 


0° samples 
USB ~~ 


=4x]1 
‘ sec (4.5) 





as well as lowering the storage requirement to: 


Satake -( 32bits | Ax 10° vas 


sample sec 
Storave S128%.10° 16x 10° 2 (4.6) 
sec sec 
MB 
Storage = 16— 
sec 


Because digital sampling was performed, there was the possibility of introducing the 
effects of aliasing. Within the software program, low-pass filtering is performed to 


eliminate the (47 f,) components, which corresponds to 2x156.8 = 313.6 MHz. This 


frequency is aliased to 1.6 MHz using the sample rate indicated in equation (4.5). As 
shown in (2.4), the output is a function of the phase difference and essentially a DC 
signal. Because of this property, the digital low-pass filter was able to utilize a very 


narrow passband. The narrow passband ensured that any aliasing introduced by the 
Do 


digital sampling would be filtered out and would not affect the final results. In theory, a 
higher decimation rate could be utilized in the implementation of this system. It would 


effectively lower the sampling rate, as well as the storage rate. 


With the decimation rate chosen, the usrp_rx_cfile.py program could be utilized 
to analyze the direction-finding system. This program captured samples and saved them 
into a *.dat formatted file. The data stored in this format could then be read and analyzed 
to determine the performance of the direction-finding system. As stated earlier, the 
phase shifter was used to implement the phase difference that would occur between a 


single signal arriving at two separate antennas. The phase shifter was capable of 


[0,180°] It was determined that 30 turns were 


implementing phase changes from 
necessary to achieve the full 180 degrees of phase change; thus, each turn of the phase 
shifter implemented six degrees of phase change. To establish a baseline for the system, 
the program was run thirty times to capture data for each turn of the phase shifter. The 
program was run for one second in order to capture an adequate amount of data and then 


saved in a *.dat formatted file. Once the data was captured, it needed to be analyzed to 


determine the differences in output caused by each turn of the phase shifter. 


The data analysis software chosen for this research was GNU Octave. GNU 
octave is a high-level language intended for numerical calculation. It provides a 
convenient command line interface using a language that is strikingly similar to 
MATLAB [4]. Octave, like GNU radio, is open-source software that is free to be used 
and improved upon by any user. It can be utilized with functions in Octave’s own 
language, which was done in this research, but can also be customized with modules 


written in C++, C, Fortran, or other languages [4]. 


The first step in establishing the baseline for the direction-finding system was to 
convert the data captured by the usrp_rx_cfile.py program into readable, analyzable data. 
This was accomplished in Octave by using the function read_complex_binary.m. This 
program opens the file specified and returns its contents as a column vector of 32 bit 
complex numbers. In research conducted before the experiment was performed, it was 


discovered that the USRP occasionally has issues with a buffer contained within the 
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FPGA. The buffer is not “flushed” as it should be when the USRP is opened or closed. 
This distorts the beginning samples of each set of data collected by usrp_rx_cfile.py. 
The open-source community is current investigating solutions to this problem, but for this 
research, a work-around was achieved. Because of the large amount of data collected in 
one second, it was possible to drop the first one million samples to remove the adverse 
affects of the non-flushed FPGA without losing too much data to make the results 
unusable. Dropping this data yielded more coherent solutions for the direction-finding 


system. 


After the data was resized, the program calculated the magnitude of the signal by 


using the formula 


Magnitude=,/I°+Q’. (4.7) 


A low-pass filter was created in Octave to filter out the high-frequency components that 
resulted from the analog mixing before entry into the USRP. After filtering and 
according to (2.8), the output of the filter should be a function of the phase difference of 
the two input signals. Within the octave program, the mean value of the filter output was 
calculated in order to determine if there was a recognizable function that resulted from 
plotting the output of the program versus the phase difference. Figure 27 shows the 


baseline established for the direction-finding system. 
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Figure 27. Relationship Between Phase Difference and Program Output. 


The experimental data is shown in circles with a best-fit curve plotted on top of 
the data. As with the hardware solution, one does not see the expected trigonometric 
relationship. While the hardware solution yielded a linear relationship, the software 
solution yields a nonlinear, parabolic relationship. It is unclear why the hardware and 
software results differ in shape, but time constraints prevented further investigation into 
this difference. One can clearly see this relationship between the phase difference and 
the program output. Using the equation for the best-fit line, the relationship between 


phase and program output is shown below: 


phase(deg) = 189. 1(, foutput —1.4 + 0.454) 
or (4.8) 
phase(deg) = -189.1( output —1.4 - 0.454) 
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Because of the above relationships, there are two solutions for every output value, which 


is a weakness of the software solution as compared with the hardware solution. 


This phase information is then entered into another Octave program which uses 
(2.4), derived in Chapter IJ, to arrive at a solution set of possible Angle of Arrivals for 
Antenna System 1. Due to the parabolic shape seen in Figure 26, each output from the 
program can lead to two possible phase angles. As shown in the previous chapter, each 
phase leads to two possible bearings. Thus, with the software solution, there are four 
possible bearings, even when using a second antenna system. This is worse than the 


hardware solution, for there is twice as much bearing ambiguity. 


Thus, the software solution has the capability of solving the problem for which 
the research was conducted. It has proved that a functioning system can be built with 
predominantly software components that can accurately determine the Angle of Arrival 
of an incoming bridge-to-bridge signal if the phase difference between a single signal 


arriving at different locations is known. 
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Vv. CONCLUSION 


A. CONCLUSIONS 
1. Performance Conclusions 


This research has shown that it is possible to perform direction finding of a 
bridge-to-bridge VHF band signal using both traditional hardware components as well as 
software-defined radio. Common techniques for conducting geolocation focused on 
Frequency Difference of Arrival and Time Difference of Arrival. This research 
confirmed that in addition to the above techniques, interferometry based on phase 
difference is a powerful tool that can also be utilized. While the software-based solution 
has the most room for future growth, the scope of this research has shown the analog 
solution to be the preferred option. Both the hardware and software solutions were 
unable to develop a single Angle of Arrival solution due to the multivalued nature of the 
arcsine and arccosine. However, the linear nature of the hardware solution and its tighter 
best-fit line yield a more reliable, preferable solution. Within the range[0,180], there is a 
unique solution for each phase. However, the parabolic shape of the software solution 
has the possibility of two phases for each output of the software system. In both these 
cases, the addition of a second antenna system greatly aids performance. Even in this 
scenario, however, the hardware solution still outperforms the software solution. For as 
stated in the previous chapter, the hardware solution yields two possible bearings for the 
source of the bridge-to-bridge signal, but the software solution yields four possible 
bearings. Thus, the bearing ambiguity in the software solution is worse than the 


hardware solution. 


There is limited room for improvement in the hardware solution. All components 
used in the experiment were top-of-the-line, standard components used in commercial 
applications. However, the hardware used in the software solution could be improved. 
While still more expensive than the hardware solution, the USRP is considered a low-cost 
device. However, it is not the only platform that can be used for software-defined radio. 


It is popular due to its affordability and flexibility. If this system were to be developed, a 
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more specialized hardware may be desired. If one is solely concerned about the bridge- 
to-bridge frequency range or other VHF signals, one would not need hardware that has as 
large a frequency range as the USRP. The USRP has capability up to the gigahertz range, 
which would not be needed in a specialized system. The specialized system could also 
utilize a much lower sample rate. Due to the slow-changing nature of the phase between 
two maritime vessels, the high sample rate is not a mandatory requirement. Although it 
would likely add cost, performance could likely be improved if hardware were utilized 


that was specially tailored for the specific, desired frequency range. 


2. Cost Conclusions 


In the design of any engineering system, the performance of the system is only 
one consideration in the decision for further development, procurement, and production. 
It is a very important factor, but in the majority of situations, the lifecycle cost is the most 
important factor. One may achieve unmatched performance, but if it is extremely 
expensive throughout the system lifecycle it will not be produced and development may 
be stopped. On the other hand, the cheapest solution may not achieve acceptable 
performance. A trade-off study must often be conducted to achieve the optimal balance 
between cost and performance. The cost associated with the hardware solution explored 


in this research breaks down as follows: 
e $220.22 for circuit board components 
e $130.00 for two circuit boards (one per antenna system) 


This does not include the test equipment utilized in the Naval Postgraduate School 
laboratory, which included the phase shifter, two signal generators, and many connecting 


wires and cables. The cost required for the software solution is as follows (from [25]): 
° $700.00 for USRP package 
° $75.00 per daughterboard 


Based on both cost and performance in the experiments conducted, the hardware solution 


is the preferable choice when compared to the software solution. However, the software 
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solution has the most room for growth. The results of future work and experimentation 
may improve the performance of the software solution to the point where it may be worth 


the extra cost. 


B. RECOMMENDATIONS 


There is significant room for future work and experimentation with both the 
software and hardware solutions. For the software solution, future work should be 
undertaken to implement more of the functionality in software. Time constraints 
precluded implementing the signal mixing in software. Additionally, a bandpass filter 
should be implemented in software and/or hardware to ensure that the signals that enter 
into the digital mixer are the frequency desired. The software experiment conducted used 
a signal generator to ensure just one frequency was analyzed by the system. However, in 
a real-world scenario, there will be many signals of many different frequencies being 
transmitted. The software would need the bandpass filter to ensure that only frequencies 
of interest are captured and analyzed. Further research could also be conducted into 
conducting the software solution experiment differently. Instead of using a decimation 
factor to lower the USB sampling rate, it could be possible to divide the information from 
the USRP into manageably-sized blocks of data. This data could then be transferred to 
the personal computer without having to sample, which would ensure that no data 


samples are excluded during analysis. 


The local oscillator addition to the hardware solution ensured that this problem 
was accounted for in the hardware solution. Thus, the hardware solution is farther along 
in the development process than the software solution. However, future work remains for 
the hardware solution. The next step in research would be to test both systems with 
wireless signals. Using hard-wired signal generators to simulate the bridge-to-bridge 
signals minimizes the loss that would be inherent in real-world signals. The noise and 
interference that would be experienced with real-world signals needs to be taken into 
account. This most likely would require the addition of a Low Noise Amplifier to both 
the hardware and software solutions. Automatic Gain Control would also most likely 


need to be implemented to ensure a consistent input voltage and power level to the 
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system. The wireless system should first be tested in a controlled environment, such as 
an anechoic chamber to shield the RF waves. After successful testing in the controlled 
environment, testing should be conducted in a real-world environment. Finally, full- 
shipboard testing should occur in a maritime environment for a final demonstration of the 


direction-finding system. 


The potential exists for additions to the system that would make it even more 
powerful to the user. Further investigation could be conducted to determine the precise 
reason why the hardware and software solutions yielded such different results. A display 
could be designed that would convert the calculated Angles of Arrival into a graphical 
display. In the current hardware system, the voltage output is supplied to a MATLAB 
program to calculate the Angle of Arrival. A device could be designed and built that 
used the voltage output and calculated the bearings without any user input. Further 
research should also be conducted to resolve the ambiguity problem. Possible solutions 
may include adding another set of antennas or perhaps comparing the two antenna signals 
in a different way that would yield only one line of bearing. To convert the results of this 
research to a final, finished product would require the testing detailed in the previous 
paragraph. It would also require the integration of the system built in this research, with 
a microprocessor that performs all calculations and a textual or graphical display for 


enhanced usability. 
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